Nanoscale ferromagnetism in non-magnetic doped semiconductors 
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While ferromagnetism at relatively high temperatures is seen in diluted magnetic semiconductors 
such as Gai-jMnxAs, doped semiconductors without magnetic ions have not shown evidence for 
ferromagnetism. Using a generalized disordered Hubbard model designed to characterize hydrogenic 
centers in semiconductors, we find that such systems may also exhibit a ferromagnetic ground state, 
at least on the nanoscale. This is found most clearly in a regime inaccessible to bulk systems, but 
attainable in quantum dots as well as heterostructures. We present numerical results demonstrating 
the occurrence of high spin ground states in both lattice and positionally disordered systems. We 
cxamine how the magnetic phases are affected by characteristics of real doped semiconductors, such 
as positional disorder and electron-hole asymmetry. 
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The Hubbard model [l| is perhaps the simplest model 
containing the essential character of correlated electrons; 
it consists of tight binding on a lattice with a purely 
on-site electron-electron repulsion. It has been studied 
extensively, e.g., in the large-t/ regime Q, on different 
latticesQ, with multiple possibly degenerate bandsQ, 
and with binary alloy disorder [fj. It has been used to 
model Mott-insulator oxides @] , high-T c superconductors 
0, organic conductors @] as well as hydrogenic centers in 
doped semiconductors [9(. For the latter, it is especially 
relevant in the insulating phase, where the Coulomb in- 
teraction is large compared to the kinetic energy. 

The Hubbard model on a lattice is defined by the 
Hamiltonian: 
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where t and U are the kinetic and Coulomb energy pa- 
rameters respectively, c\ a {ci a ) is the usuai electron cre- 
ation (annihilation) operator on site i with spin a, and 
angular brackets denote nearest neighbors. For example, 
in the tight-binding model with hydrogenic wavefunc- 
tions, t(r) = 2(1 + r/ae) exp(—r/aB)[lO(. Nagaoka 11 
showed that in the limit U/t — > oo, the Hubbard model 
on a finite bipartite lattice of dimension d > 2 with pe- 
riodic boundary conditions and a single hole (away from 
half filling), has a ferromagnetic ground state. The rea- 
son for this can be understood by considcring a system 
with hole density 5. The kinetic energy gain of the holes 
is more restricted in an antiferromagnetic background of 
spins than in a fully spin-polarized background [12J. This 
restriction is of order St, whereas the antiferromagnetic 
superexchange cost is order J = 4i 2 /[/[13]. At small 
S and large enough U, 5t > J, and the system prefers a 
ferromagnetic configuration because it allows for less con- 
fined carriere. Subsequent work to Nagaoka's pioneering 
result has shown that ferromagnetism is a subtle effect 
depending on lattice geometry. For example, Lieb and 
MattisjlJ] proved that in finite one-dimensional systems 
with zero-wavefunction or zero-derivative boundary con- 



ditions, the gro und state must be a singlet, and Haerter 
and Shastry[15( recently showed that on the frustrated 
triangular lattice (t < 0) , an itinerant hole actually hclps 
to produce an antiferromagnetic ground state. 




FIG. 1: Zero temperature mean-field theory phase diagram 
of the Hubbard model on a simple cubie lattice (512 sites). 



The mean field diagram of this model on a simple cubie 
lattice is given in figure[lj which agrees qualitatively with 
the more extensive work in two dimensions by Hirsch[16j. 
An antiferromagnetic phase exists at half-filling for ali 
values of U/t, due to the the effective antiferromagnetic 
interaction from superexchange. Such a mean-field anal- 
ysis does not include the possibility of phase separation, 
e.g. the existence of polarons corresponding to "carrier- 
rich" ferromagnetic and "carrier-poor" antiferromagnetic 
spatial regions. Phase separation may alter the simple 
phase diagram of fig. [TI substantially [1 7j . though Dagotto 
et al. [IH argue for its absence in the Hubbard model 
based on 10- and 16-site square lattices. Even though 
its precise location in phase space depends on dimen- 
sion and requires more careful work, it is clcar that for 
large enough U/twe expect an antiferromagnetic to ferro- 
magnetic transition on some mesoscopic or macroscopic 
length scale, as a function of doping. 

However, in spite this expectation, Nagaoka ferromag- 
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netism does not seem to have been seen experimentally. 
Whereas most Mott-insulator oxides appear not to have 
a large enough U/t ratio to allow for it, doped semi- 
conductors have U jt tunable over several orders of mag- 
nitude (due to the exponential dependence of the hop- 
ping t on the dopant spacing), and thus are a promising 
candidate for Nagaoka ferromagnetism. Unfortunately, 
the dopants are not arranged on a superlattice, and the 
Hubbard Hamiltonian must be modified to include a site- 
dependent hopping term. 
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It has been shown[l|| that in uncompensated semicon- 
ductors (half-filled Hubbard band), the randomness of 
the dopants results in a valence-bond glass (random sin- 
glet) phase being a better description of the ground state 
than the antiferromagnet predicted on a lattice. Exper- 
iments on compensated systems (away from half-filling), 
show no evidence of ferromagnetism in conventional semi- 
conductors with non-magnetic dopants [2(J. This can be 
attributed to the localization of the holes due to the 
strong randomness in hopping parameters that results 
from random dopant positions. The holes are conse- 
quently unable to gain the kinetic energy which favors 
a spin-polarized background. Thus, despite the ability 
to tune U/t over such a large range, Nagaoka ferromag- 
netism remains elusive. 

In this work, we show that there exists a regime 
of doped semiconductor systems that is attainable in 
nanoscale quantum dots and heterostructures, but not 
accessible to bulk systems, which is more suited to the 
occurrence of Nagaoka type ferromagnetism. In this re- 
gion, ferromagnetism is found at least at the nanoscale, 
and has a higher likelihood of emerging on macroscopic 
scales (e.g. in modulation doped systems). 

The existence of such a regime is suggested by special 
properties of the hydrogen atom (and hydrogenic centers) 
that effectively reduce the disorder and move the ran- 
domly doped semiconductor in the direction of a lattice 
problem. The key property is that the second electron 
of a H~ ion is bound by only 0.0555 Ryd[2l|, an energy 
much less than the 1 Ryd binding of the initial electron. 
This property is linked with the fact that the two-electron 
wavefunction of H ~ is spatially much larger than the ls 
wavefunction of the hydrogen atom. It is much easier for 
the second electron on a hydrogenic center to hop away 
than it is for a single electron on such a center to do so. 
Thus, near half filling, in Hubbard model parlance, the 
hopping amplitude for an electron is much larger than 
for a hole. Consequently, starting from the half-fillcd 
system (i.e., the uncompensated doped semiconductor), 
the system with a small percentage of extra electrons, 
because of the more extended wavefunction, experiences 
a greatly reduced effect of the positional disorder as com- 
pared with the corresponding hole-doped (less than half 



filled, i. e. compensated) system, and so will behave more 
like the uniform lattice. In diluted magnetic semiconduc- 
tors (DMS), the existence of the relatively larger Bohr 
radius of the carriers (wlOÀ) compared to the extent of 
the spins on the magnetic ions (1-2À) allows the carrier- 
magnetic moment interaction to dominate, resulting in 
a ferromagnetic ground state [22]. In the clcctron-dopcd 
semiconductor, electrons occupying the D~ state have a 
larger Bohr radius than the electrons giving rise to the 
effective exchange interaction (J ~ t 2 /U). This could 
cause carrier hopping to dominate, similarly resulting in 
a ferromagnetic ground state. At the very least, the dif- 
ferent radii of the doubly vs. singly occupied sites sug- 
gests that we modify the Hamiltonian ((2]) to become: 
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(tij(rii,n j )cl T c j(7 + ce) + U^n l] n ii (3) 



where rij is the total occupation of site i, and iy is now 
has an occupation dependence given by: 



tij (jli , Tlj ) — 



t rij ■ — 2 and n% = 1 
t otherwise 



where t is larger (and can be much larger) than t. For 
greater than one electron per site, the low energy spec- 
trum in the limit U 3> t, t is given by the t — J model, 
where J = 4t 2 /U, as one would expect 23j. Hirsch has 
investigated a similar Hubbard model with occupation- 
dependent hopping, but with focus on its prediction of 
superconductivity pairing 2J| . We proceed with semicon- 
ductors in mind, motivated by the notion that one might 
expect to find a ferromagnetic phase in this model. 

In the present work, we focus on the effect of changing 
the ratio t/t on small random clusters and lattice sys- 
tems. We have computed U/t and t/t appropriate for 
dopants in semiconductors by performing a rcalistic cal- 
culation of single particle states of donors placed on a 
simple cubie lattice. We extract the dependence of these 
ratios on the simple cubie lattice Constant by fitting the 
single particle bands to a tight binding model. Details of 
this work will be given elsewhere 25]; here we show the 
resulting parameter ratios in figure [2] The Mott metal- 
insulator transition criterion in these units is R/a^ = 4. 
We see clearly that the range of U/t and t/t can be var- 
ied substantially in doped semiconductors. The large 
span of U /t originates in the exponential dependence of 
the hopping parameter on the atomic spacing, and the 
variation of t/t from the relatively large size of the two- 
electron wavefunction appearing as a factor in this expo- 
nential. Here we restrict ourselves to ratios in the ranges: 
U/t = [10, 100] and t/t = [1, 10], which are conservative 
when compared to the physically attainable ranges. 
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FIG. 2: Variation of ratios U/t and t/t with the dopant spac- 
ing (related to the dopant density p by p — -gj. 

We have solved the Hubbard and corresponding t-J 
models exactly using the Lanczos method to diagonalize 
the Hamiltonian matrix, after first taking advantage of 
ali available spatial symmetries of the system and the 
spin symmetries (where allowed by memory constraints). 
The Hubbard model depends on U/t and t/t, whereas the 
t — J model only depends on t/J = j(t/t)(U/t). Thus, 
the value of t/J marking the onset of the Nagaoka state 
defines a line in U/t vs. t/t space. Figure [3] shows the 
ground state spin phase diagram for the 8-, 10-, and 16-site 
square lattices doped with one extra electron. One sees 
that the increase in t/t causes the area of the maximal 
spin ground state phase to increase for electron-doped 
systems. 




t/t 

FIG. 3: Ground state spin diagram resulting from the exact 
diagonalization of © on 8-, 10-, and 16-site square lattices 
with 9, 11, and 17 electrons respectively. Hubbard model re- 
sults are displayed as open symbols. Closed symbols and lines 
show the result of the corresponding t — J model as described 
in the text. S ma x is the largest allowed spin on each lattice, 
and Siow denotes a low ground state spin (1/2 for 10- and 
16-site clusters, the minimum value). 



We have also solved the Hubbard model on square lat- 
tices with a single hole, and see virtually no dependence 
of the ground state spin on t/t, indicating a pronounccd 
electron-hole asymmetry. It is clear that this asymmc- 
try originates from the electronic states having greater 
radius than the hole states, since the square lattice is bi- 
partite and thus for t = t the problem is electron-hole 
symmetric. In the case of a non-bipartite lattice or clus- 
tcr, there is an asymmetry even in the simple Hubbard 
model[26j], and taking t > t may yield a combined effect 
whereby the initial asymmetry is enhanced. 

Figure H] shows several clusters, each less symmetric 
than a lattice but retaining some spatial symmetries, 
that have a high-spin ground state in certain parameter 
ranges. The magnetic phase diagram of the two largest 
clusters are shown in figure OH where t/t is plotted vs. 
a parameter measuring the geometry of the cluster. We 
again see that the increase in t/t enlarges the high-spin 
rcgion of the phase diagram and thus indicates that the 
high-spin state becomes more robust. 
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FIG. 4: Geometries of several high-spin clusters. The spin 
given is the ground state spin when U/t\ = 100 and £2/^1 = 
0.5. These high spin ground states persist in considerable 
regions of phase space around this point. 



To further test the robustness of a cluster's high spin 
ground state, we multiply each ti by a random factor 
A whose logarithm is chosen from the box distribution 
P(logA) = l/(21oga), logA e [-Ioga, + log a). This 
more accurately characterizes the fìuctuations we expect 
in a random system, since the hopping is exponentially 
dependent on the inter-site distance and we do not expect 
the fìuctuations to preserve any symmetry present in the 
cluster. We start with a cluster known to have a large 
spin ground state (e. g any cluster in fig. HJ) and average 
over many of the just described random perturbations. 
We find that the percentage of randomly perturbed clus- 
ters that retain the high spin ground state of the originai 
clustcr increases dramatically with t/t (set to the same 
value on ali edges). Table U shows how this percentage 
depends on t/t and a for the particular cluster IV of 
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FIG. 5: Ground state spin diagram for cluster IV (left) and 
V (tight) of fig.[4]with 10e~ as a function of t/t. Lines denote 
U = 100 (solid), U = 50 (dashed), and U = 20 (dash-dot). 



fig.Hwith ti = 1, i 2 = 0.3, J7 = 100. We see clearly that 
increasing t/t makes the high spin ground states of the 
clustcrs in fig. [5] significantly more robust to geometrie 
fluctuations present in the actual system. 
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TABLE I: Percentage of random instances of cluster IV (fig. [4} 
that retain a high spin (S = 3) ground state. Clusters are 
obtained by starting with t\ = 1, Ì2 = 0.3, U = 100, and mul- 
tiplying t and t on each edge by a random factor in [a,l/a]. 



In conclusion, we have shown that on both finite (pe- 
riodic) lattices and several less symmetric clusters, in- 
creasing t/t makes the appearance of nanoscale ferro- 
magnetism significantly more likely. The effect is greater 
for large U/t, which implies small dopant densities, well 
below the metal-insulator transition density. Such high 
spin states should be observable in doped quantum dots 
with dopant number Nd = 6 — 15 and with a small excess 
of electrons N e — Nd = 1 — 2. The same density and 
(excess) electron doping regime is also the most likely re- 
gion for the possible appearance of true macroscopic fer- 
romagnetism,e.<?. in modulated structures with dopants 
in both quantum wells and barrier regions, so that the 
region of excess electrons can be achieved, unlike in a 
true bulk doped semiconductor. However, obtaining a 
conclusive answer to this question numerically requires 
going beyond the small sizes possible with exact diago- 
nalization methods, using, e.g., density matrix or pertur- 
bative renormalization group methods, in combination 
with other numerical techniques. Even if true ferromag- 
netism on the macroscopic scale is absent, our calcula- 
tions show that there should be a significant asymmetry 
between the magnetic response of systems with excess 
electrons above the half-filled (uncompensated) case, and 



those with a deficit of electrons from the half-filled case, 
(Le. traditional compensated) : the former should have 
a larger susceptibility in the paramagnetic phase at low 
temperatures. This can be experimentally checked by us- 
ing gates to tune the electron density. If ferromagnetism 
is attained on large enough length scales, it may show up 
as hysteresis in transport measurements due to magnetic 
domains. 

This research was supported by NSF-MRSEC, Grant 
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